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I. INTRODUCTION 


The possible existence of strange quark stars which are made entirely of u, d and s quarks is 

It has been conjectured that the 


B- 


one of the most exciting aspects of modern astrophysics Q|-P|. 

compact objects such as X-ray pulsar Her X-1, X-ray burstar 4U1820-30 are likely strange stars[ 
P], In general, perhaps pulsars can be divided into two groups: one is modeled as neutron stars, 
and the other as strange stars. It is essential to distinguish these two candidates of pulsars from 
observations and to predict their properties furthermore. 

Usually, if strange stars do exist, the basic differences between strange stars and neutron stars 
are their dynamical properties and thermodynamical behaviors. For example, firs^, the mass- 
radius (M-R) relation of strange stars is quite different from that of neutron starsespecially, 
this relation depends on the rotation frequency of compact object. Secondly, for rotating stars, 
the moment of inertia plays an important role to investigate their dynamical and electromagnetic 
behavior. The frequency of radio signals emitted from pulsars contain much information on the 
moment of inertia of the sources, which helps us not only to identify the type of pulsars - neutron 
stars or strange stars, but also to understand some phenomena and mechanism. It is generally 
accepted that the sudden spin-ups, glitches of pulsars, can be explained from their angular mo¬ 
mentum transfer between their crusts and inner fluid. Obviously, this transfer depends on their 


moments of inertia 


Thirdly, a newly born neutron star or strange star may pass through 


various stages of early evolution. These stages are determined by neutrino time scale and this 


scale is proportional to the star radius and the neutrino mean free path [3]. These two physical 
quantities are dominated by temperature and the equation of state (EOS) remarkably. In fact, 
due to the considerable difference between the EOS of the strange quark matter and the neutron 
matter, we can distinguish the strange star and the neutron star from observation. The EOS at 
finite temperature plays the key role to understand the behavior and the evolution of compact 
objects. 

The early calculations of EOS for strange star are based on the MIT bag model j^. LatCT, 
many other models such as vector interaction and de nsity -dependent scalar potential model 
the quark mass density-dependent (QMDD) model 111 and etc. were employed. 


The QMDD model was first suggested by Eowler, Rata and Weiner [1^ and then used by many 
authors to study the stability and the thermodynamical properties of strange quark matter |l^ - 
The basic hypothesis of the QMDD model is that the masses of u, d quarks and strange 
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quarks (and the corresponding anti-quarks) are given by 


mq = -— {q = u,d,u,d), 

Sub 


(1) 


ms,s = ruso + 


B 

Sub' 


( 2 ) 


where ns is the baryon number density, niso is the current mass of the strange quark and B is 
the vacuum energy density. Eqs. m and m can easily be understood from the quark confinement 
mechanism, if we notice that the density approaches to zero and the masses of quarks tend to 
infinite when the volume of the system goes infinite . This confinement mechanism is almost the 
same as that of the MIT bag model. It has been shown by many authors that the thermodynamical 
properties given by QMDD model are similar to that of MIT bag model and the M-R relation, 
viscosity radial oscillation of both rotating and non-rotating strange quark star described by QMDD 
model are qualitatively similar to those obtained with MIT bag model Q. 

Although the QMDD model can provide a dynamical description of confinement and explain 
many aspects of strange quark matter, it suffers from a basic drawback that it cannot reproduce 
a correct lattice QCD deconhnement phase diagram because the quark masses are divergent when 
—> 0. To excite an infinite weight particle, one must pay the price for infinite energy, i.e. 
infinite temperature. It means that the QMDD model is a permanent quark confinement model. 
It is unable to describe the deconhnement process from neutron matter to quark matter inside the 
pulsar. To overcome this difficulty, in a series of previous papers Q-Q we suggested a quark mass 
density- and temperature- dependent (QMDTD) model which is based on a non-permanent quark 
conhnement Friedberg-Lee model and proved that not only the above difficulty can be overcome 
but also many physical prop erties, for example, the stability, the deconhnement phase transition 
of strange quark matter 17|- 1^. the binding energy of dibaryon system can be explained Q. 
Employing QMDTD model, Gupta and his co-workers studied the M-R relation and the radial 
oscillations of proto strange stars and found a lot of reasonable and interesting results But 
their study was limited to non-rotating proto strange stars only. 

This paper evolves from an attempt to extend the study of Gupta et. al. to slowly rotating 
proto strange stars by using the QMDTD model. We will employ the Hartle’s method j22|| to sketch 
the main features about strange stars rotating with the frequency much smaller than the Keplerian 
limit. We will discuss two cases: (1) Uniform rotation where angular velocity D = const.; (2) D 
is a function of r, D = D(r), which will be called differential rotation below. We will investigate 
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the temperature dependence of M-R relation, of the moment of inertia and of the effect of frame 
dragging for proto rotating strange star. 

The organization of this paper is as follows. In the following section, we will give a brief review 
of the QMDTD model. By means of the Hartle’s formalism, we will discuss the slowly rotating 
strange star with uniform angular velocity in Sec.Ill, and that with differential rotation in Sec.IV, 
respectively. The last section contains a summary. 


II. THE QUARK MASS DENSITY- AND TEMPERATURE- DEPENDENT MODEL 


rrief review of the QMDTD model in this section. The detail of this model can be 

1 

19l| . Here we only write down the main steps which are necessary for calculating 
the EOS for strange quark matter. 

According to the QMDTD model, the masses of u, d and s quarks (and the corresponding 
anti-quarks) are given by 


We will give a 

n 

found in refs. 171 - 


rug 


-- {q = u,d,u,d), 

SriB 


( 3 ) 


rus^s = rnso + 


B{T) ^ 

Sub 


( 4 ) 


Comparing Eqs. and 0 with Eqs. o and 0, the constant B in QMDD model is replaced by 


B{T) = Bo[l-{T/T,f], 0<T<T, 


( 5 ) 


B{T) = 0, T>n (6) 

where Bq is the bag constant at zero temperature and Tc = 170MeU is the critical temperature 
of quark deconfinement phase transition. The basic extension of QMDTD model is that B{T) 
depends on temperature because according to the Eriedberg-Lee soliton bag model, the vacuum 
energy density of the bag equals the different value between the local false vacuum minimum and 
the absolute real vacuum minimum and this value depends on temperature Q 
The thermodynamic potential reads 

Z = -S"t[ ln(l + (7) 

V Jo dk 
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Figure 1: Diagram for the EOS of the strange matter in QMDTD model. 


where i stands for u, d, s (or u, d, s) quarks, is the corresponding chemical potential (for anti¬ 
quark = —jJi), £i = y^rn?”+^ is the single particle energy and is the mass for quarks and 
anti-quarks. dNi{k)/dk is the density of states for various flavour quarks. The density of states of 
a spherical star has been calculated in ref.j^ 


Ni{k) = Ai{kRf + Bi{kRf + Ci{kR), 


( 8 ) 
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( 9 ) 


A- = — 

* Ovr ’ 


^ % rr. 1/ N /"llN TTt 


^ (7j 1 , /c mj, 1, A: , ttA: 

= |^{^ + — + i:^ tan-^— - — 
Ztt 6 rrij k rrij Irrij 


( 10 ) 


+( 


m- 


M1.45 


9i 




( 11 ) 


k' 3 _ 42 (!^ _ 6 . 5)2 + 100 - 
where gi is the total degeneracy. For a slowly rotating star, as a zero-th order approximation of 
the Hartle’s perturbation formalism, we use the density of states of a spherical cavity to calculate. 

After getting the therm ody namic potential Z, we can derive the energy density p and the 
pressure p. The results are 1^ T^- 20|) 


p = 


1 d{Q/nB) 


V d{l/nB) 


n riR cin 

= -77 + 




V V duB 


( 12 ) 


Tyfli 


n ^ T fdn\ 

P=-i7 + T".«-7(^j 


( 13 ) 




The extra terms in Eqs. m and (d come from the dependence of the quark mass on the baryon 
density. 

The EOS for strange quark star can be calculated from above formula and the result is shown in 
FiglU where the parameters are fixed as Bq = IJOMeVfm niso = 150MeF and = ITOMeF. 

III. UNIFORM ROTATING STRANGE STARS 

A. Formalism 


In this section, we employ Hartle’s formalism to study the slowly rotating strange star. We 
review this method briefly below and the details can be found in ref. [2^. 

The metric of non-rotating configuration is 

ds^ = + r'^idO^ + sin^ edp^). 

In hydrostatic equilibrium, i^(r) and A(r) are determined by the TOV equations 

2 dP 


dv 

dr 


p + P dr ' 


(14) 


(15) 
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(16) 


gAM = (1 _ ^) 


dm , 9 

—— = 47rr p, 
dr 


(17) 


dP 

dr 


{P + p){m + Anr^P) 


(18) 


r2(l-^) 

where P and p are the pressure and energy density of strange matter respectively and they relate 
each other by EOS. The stellar mass M = m{R), where R is the radius of the star and is obtained 
by solving equation P{R) = 0. For a slowly rotating star, the configuration is no longer static and 
spherically symmetric, but axially symmetric instead. The metric becomes 




2m 


-]dr^ 


r — 2m 

+r^(l + 2k)[d9‘^ + sin^ 9(dip — Lodt)'^] + 0{io^), 


(19) 


where uj{r) is the angular velocity of a locally non-rotating inertial frame. Expanding the metric 
in spherical harmonics, we find 

h = h{r,9) = ho{r) + hi{r)Pi{9) + h2P{9) + ... 
m = m(r, 9) = mQ{r) + mi(r)Pi(0) -|- m2P{9) + ... (20) 

k = k{r, 9) = kQ{r) + ki{r)Pi{9) + k2P{9) + ... 

It can be proved from symmetric consideration that hi{r), mi(r) and ki{r) vanish. Therefore, 
under the slow rotation approximation the metric can be written as 


ds^ — —e^[l + 2(/io + h2P2y\dt^ + e^[l + 


2(mo + m 2 P 2 ) 


r — 2m 

-|-r^(l + 2k2P2)[d9'^ + sin^ 9{dt — ujdip)^] + 0{uj^) 


]dr‘^ 


( 21 ) 


where P' and can be calculated by TOV equations (|T5)) - ((T8)) and P 2 = {3cos9 — l)/2. 

The rotation not only affects the metric, but also can change the distribution of pressure in the 
interior of the star. In a co-moving reference frame with the fluid, the pressure correction is 


and the energy density correction 


P+{P + p)p* = P + AP, 


P + {P + p)^p* = p + Ap. 


( 22 ) 


(23) 
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Here p* indicates the dimensionless corrections on the pressure. It can also be expanded as 


p*{r,e) =Po{r) +p*2{r)P2 + ... 


The stress-energy tensor in a rotating star is 

tWi' + AP)u^u‘' + (P + AP)j;;, 


(24) 


(25) 


where 


= {-gtt - 


= VLu\vr =u^ = 0 


(26) 


H is the angular velocity of fluid observed far from the star. Noting that the perturbed Einstein 
equation + SttAT/^ and the zeroth order equations we get 


ag;: = sttAt;. 


For the component of t = and /r = (/> of Eq. (EH) , we have 


I d . dw Adj e d 


dr 

1 d , d 

= ~7rv J Gtt} + 0 . 3 . 7 ^ 

dr dr r^ sin^ 0 dd 


r dr 

dn 

dr 


<™ "tm*- 


where 


j{r) 


— g ~('^+'^)/2 


(27) 


(28) 


(29) 


and we have defined w = H —a;. For the uniform rotation, the angular velocity of fluid is a constant 
Q = const.. Therefore, the right hand side of Eq. (EHI) vanishes. We expand tz7(r, 6») Q: 

( O') ^ 

^{r,0) = > tz7;(r)(— 

s 

1=1 

Then the radial function of wi{r) satisfies 


sin0 d9 


(oiA, + df - -b n., = 0, 

dr dr r dr 


At large r, 'cui{r) has the form 


(30) 


(31) 


UJi{r) const.r ^ ^ + const.r 


J-i 


(32) 







At infinity where the spacetime is flat, lo must decrease faster than 1/r^. Because oi w = Q. — u, 
only the coefficient of / = 1 in the w expansion does not vanish and Eq. (EU reduces to 


I d ^ dw A dj 

— — {r — } + -—w = 0, 
dr dr r dr 

where w = TUi{r). Outside the star j(r) = 1, the solution reads 


(33) 


(34) 


where the constant J is identified as the total angular momentum of star and satisfies 

J = I r^^/^gdV, (35) 

where is the determinant of metric of space components. The moment of inertia is simply 

I = ( 36 ) 

The components u = g, = t and p = g = r ol Eq. (EH) and the terms / = 0 in the expansion (Oni) 
lead to 


d'diQ A_ 2 */ , TD\^P \ ^ -2 4 f dio\ 1 3 ^ 7 ^ 2 


dr 


= 47rr p^{p + P)— + —j r 


r -j—'ca , 


dP 12 \dr J 3 dr 


(37) 


dho , 1 ^ A7r{p + P)r‘^ ^ 1 .2 f dujy 

P^^P+—) = — -TZ—1^0-TZI— —d ( — ) • (38) 


dr (r — 2m)^ ' r‘^' r — 2m \2r — 2m'' \dr J 

Another constraint between the perturbations of the metric and those of energy and pressure 
comes from the hydrostatic equilibrium, i.e. the conservation of stress-energy tensor of perfect 
fluid Py = 0. In the non-rotating case, it can be expressed as 


Pi 12 i 

- ^ = 0 


p + p 2v 

and in the rotating case, the corresponding perturbation equation reads 

1 


P*i + h,i - sin^ 0),* = 0. 


(39) 


(40) 


For the terms / = 0, we find 


dr dr 3 dr 


Combining Eqs. m and (TTHl) . we find 


dPo 


-f 


1 


dr ' 12 r — 2m 
A'k{p -|- P)r‘^ 


J 


f dijj\‘^ 


1 d ,r^rvj‘^. 


-(- 


r — 2m 


Po + 


3 dr r — 2m 

fhnr'^ , „ 1 , 

“ (SttPT^). 


(r — 2m)^ 


(41) 


(42) 
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Figure 2: The radius R in km vs. mass M/Mg at the temperature T = OMeV for different angular velocity 

n. 


On the other hand, one can prove that the increasing of stellar mass 5M and mean radius 6R are 
given by 


6M = mo{R) + 


(43) 


6R = pI{P + p)/ . 
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(44) 










Figure 3: The radius R in km vs. mass M/Mg at the temperature T = bOMeV for different angular velocity 

n. 

When EOS, 0 and the boundary value of P{0) are given, we can integrate the TOV equations 
until reaching the surface of the star P{R) = 0. Then we solve the Eq.(jSSl), Eq.(jSZl) and Eo. dl^ 
to get the perturbations and finally calculate the moment of inertia by Eq. p6|) . the correction of 
the mass and radius in rotational case by Eq. (01 and Eq. 01 ). We do this numerical calculation 
with the rotational frequencies smaller than the Keplerian frequency 

0^ < _ M/R^, (45) 
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Figure 4: The increase of maximum mass 5Mmax/Ms vs. the angular velocity O at different temperatures. 


which is 


required by the Hartle’s method 
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B. Results 


1. M-R curves 

The Mass-Radius relation is one of the main results which distinguishes the strange stars and 
the neutron stars. In our calculation, the M-R relations not only depend on temperature T, but 
also on angular velocity fl. Fix the temperature T = 0 and 50MeV, the M-R curves for different 
Q are shown in FigEl and 0 respectively. By using FigEl FigEland Ea. MfiB . we find Keplerian 
frequency ~ 1.5 x lO'^Hz, which is much larger than the frequencies used in FigEl and FigEl 
It is seen that the maximum mass and the corresponding radius increase with Q. To show the 
temperature effect more transparently, we plot the M^ax vs. curves for different temperatures 
in FigEl and see that for a fixed fl, Mmax increases with temperature. The effect of rotation 
becomes more and more important when temperature increases. From FigEl we can express the 
increase of the maximum mass by a simple formula 

OC 0^, (46) 

where Mg is the mass of the sun. This is very reasonable because all perturbations in Hartle’s 
method are assumed to be proportional to the square of the angular velocity of the fluid. 

2. The effect of frame dragging 

The frame dragging is a general relativistic effect for a rotating object. To show the effect 
of frame dragging in a slowly rotating star, we plot the uj{r)/fl vs. r/R curves for different 
temperatures in FigEl We see from FigEl that this effect decreases when one approaches to the 
surface of star, and the dragging increases when temperature increases. 

3. The moment of inertia 

Besides the M-R relation, much information for the configuration of a strange star is provided 
by its moment of inertia, which is one of the essential facts to understand many observational 
phenomenons. FigEl shows the relation of the moment of inertia-Mass. We see that the moment 
of inertia will increase with temperature. 
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Figure 5: The ratio of angular velocity of inertial frame and fluid in the interior of the star uj{r)/VL vs. r/R, 
where R is the radius of the star, at the temperatures T = 20, 5QMeV. 


4- Horizontal-branch oscillation 


Recently Stella and Vietri 


2i\ made a link between quasi-periodic oscillations(QPO) and general 


relativistic Lense-Thirring precession caused by frame dragging. The low frequency of QPO, called 
horizontal-branch oscillations(HBO), is the frequency of the L-T precession of an inclined circular 
orbit with Keplerian orbital frequency. Psaltis et al. ^ found that the ratios of the moments of 
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Figure 6: The moment of inertia 745 vs. mass M/Mg at the different temperatures T = 0, 20,50, iQMeV for 
the star rotating at the frequency = 200i7z. 

inertia and masses 145 ji-^) of five bright sources (GX 17+2, Cyg X-2, GX 5-1, Sco X-1, GX 
340+0) satisfy I^^KM/Mg) > 2.3. Kalogera and Psaltisj^ calculate this ratio for various models 
of neutron star and argue that this requirement cannot be satisfied for neutron stars. 

Here we suppose these hve sources are strange stars and discuss whether the condition about 
I 4 ^/{M/Ms) is satisfied. Fig|71tells us that to satisfy the condition I 4 ^/{M/Mg) > 2.3 one require 
too high temperature which is over 120MeV from the Figure. This makes the validity of the 
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Figure 7: The ratio of the moment of inertia and the stellar mass /Mg) vs. mass M/Mg at the 

different temperatures T = 0, 20, 50, SOMeV. 

suggestion by Stella and Vietri to be challenged also for the strange stars. 


IV. DIFFERENTIAL ROTATION STRANGE STARS 

We will study the configuration of differential rotating star, i.e. the angular velocity of the fluid 
0 = Q{r,9), in this section. Obviously, differential rotation will affect the evolution of a hot and 
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new-born star. Noting that Hartle’s method can only be used to study the uniform rotational star, 
an extension is needed when star is newly born and differentially rotating. 


A. Formulism 


It can easily be seen that the Eq. (EHI) is still valid in the case of differential rotation, but the right 
hand side does not varnish now. Expanding n(r, 0) (and then w{r,9)) by the Legendre functions 


Ea. d^Hl) becomes 


^{r,9) = ^n;(r)(- 


1 dPi 


1=1 


sin0 dO 


)• 


(47) 


dr dr r dr r^ 

P dr dr r^ 


(48) 


Remember that outside the star, there is no fluid, 14/(r) = 0, therefore, ro/(r) in Eq. (m still has 
the form 


ru/(r) —> const.r ^ + const 


J-i 


A similar discussion can be made as that of the uniform rotation case: u should decrease faster 
than 1/r^ if the spacetime is asymptotically flat. As a result, only the term with I = 1 remains, 
while ro = ■a7i(r) and 14 = 14i(r) are the functions of radial coordinate only. Eq. ( 1151 ) is simplified 
to be 


1 d ^ dw Adj 

[r j{r) — } + -—w = 


dr 


dr 


r dr 


1 d f 4^ \ dl4. 

j(r) —}, 


dr 


dr 


with w'{t)) = 14'(0). The solution is 


w{r) = const. 


^3 


(49) 


(50) 


Here the constant J can still be interpreted as the total angular momentum of the star. We obtain 


w{R) = Q{R) - ^ (51) 

at the surface of star r = R. In differential rotation, the angular velocity of fluid depends on the 
radial coordinate. This means that in the internal regime of the star, a very thin layer r ^ r + dr 
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of the fluid rotates with the same angular velocity ^{r) and generates angular momentum of this 
layer dJ, which satisfies, 


dJ = Q{r)dl. 


dl is the moment of inertia of this thin layer. Hence the total moment of inertia of the star is 


I = 


dl = j dJ/n = j 


(52) 


where the Eq. ( 1551 ) is used at the last equality. Since the differential rotation changes the derivatives 
of stress-energy tensor in hydrostatic equilibrium, we find 


P,i 


p + p w 


^ = 0 , 


(53) 


where 


9t(j> 

9tt ‘^^9tif> “ 1 “ ^‘^9<j><j> 

In the slowly rotating approximation, the effect of differential rotation is considered as perturbation. 
Therefore we have 




P*i +h^i —siTT 6e '^),i + (e sim 0)n^i = 0. 


(55) 


For the term I = 0 and i = r, Eq. (|P^ gives 

Po,r + ho,r - = 0. 

Combining with Eq. we find 

^ 47r(p + p) ^2 ^ mo(87rpr^ + 1) du ^ 

dr r — 2m (r — 2m) 12 (r — 2 m) dr 

1 d 2 r^j^w dQ 

- ( — 1 -1 J- 2 - 

3 dr V — 2m ^ 3 r — 2m dr ’ 


(56) 


(57) 


with boundary conditions mo = Po = 0 - Eo. iUTtI) is valid in the differential rotation, because the 
derivatives of that Q(r) do not appear in both of the right and left hand sides of En. (l27l) as in the 
uniform rotation. We can solve En. (1371) (j57|l uumeri cal]v. 


B. Results 


We employ the distribution of the angular velocity of the fluid 

n{r) = 


(58) 
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Figure 8: The radius R vs. mass M/Mg at the temperature T = OMeV for the different values of /3. 

where Qc is the angular velocity of the core and P is a constant, to investigate the influences of 
differential rotation. Fixing the surface angular velocity = WSAdHz, which is still smaller 
than the Mass-shedding limit at the surface of a uniformly rotating star and satisfies the slowly 
rotating approximation, we show the M-R curves for at temperature T = 0 and 50MeV in 
FiglHl and FigEl respectively. We see that the effect of /3 on Mmax is remarkable, and the result 
reduces to the uniform rotation when /3 = 0. A differentially rotating strange star can sustain more 
mass than that of the uniform rotating configuration. The curves of the maximum mass increase 
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Figure 9: The radius R vs. mass M/Mg at the temperature T = 5QMeV for the different values of /3. 

5Mmax/Mmax VS. Vtc/^e, the ratio of the angular velocity at the core and at the surface, are shown 
in FigfTUl We see that dM^ax/Mmax changes considerably with temperature and Oc/f^e- A hot 
and differentially rotating strange star will suffer more maximum mass increase. 

Similarly, we can address the effect of on the moment of inertia at temperature T = 0 
and 50MeV. We show this effect in Fig il 1 l and Figrespectively. Comparing to FiglHland FigEl 
we find that corresponding to the mass, the effect of (5 on the moment of inertia is small. This 
is reasonable and expected, because in the slowly rotating approximation, the effect of rotation is 
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Figure 10: The ratio of the maximum mass increase 5Mmax/Mmax vs. Oc/Og, the ratio of the angular 
velocity at the core and at the surface at the different temperatures T = 0, 20, 50, SOMeV. 

too weak to change the star’s shape (almost spherical), which has strong influence on the moment 
of inertia. 
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M/M 
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Figure 11: The moment of inertia 145 vs. mass M/Mg at the temperature T = OMeV for the different values 
of (3. 


V. CONCLUSION AND DISCUSSION 


As shown by Gupta et. al. 


22 l |. the parameters of the EOS in the QMDTD model have few 


influences on the configuration of the strange star. We extend their work to the rotational case 
and have studied the slowly rotating strange star with uniform angular velocity, employing the 
QMDTD model and the Hartle’s method. The M-R relation , the moment of inertia and the frame 
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Figure 12: he moment of inertia 745 vs. mass M/Mg at the temperature T = 50MeV for the different values 
of (3. 

dragging of strange star with different rotating frequencies are given. We found that mass, the 
moment of inertia and the frame dragging increase when temperature increases. We also found 
that the challenges of Stella and Vietri for HBO and the moment of inertia > 2.3) for 

five bright sources being neutron stars also exist for the strange star. We cannot use the strange 
star to solve this difficulty. 

Furthermore, we have extended the Hartle’s method to study the differential rotating strange 
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star with Q = i7(r). We found that the massive strange star can be prepared in the differential 
rotation. In the slowly rotating approximation, compared to the temperature the differential 
rotation has less effect on changing the moment of inertia of the star. 
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